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Abstract 

We consider a generalized two-species population dynamic model and analytically solve it for the amensalism and 
commensalism ecological interactions. These two-species models can be simplified to a one-species model with a time 
dependent extrinsic growth factor. With a one-species model with an effective carrying capacity one is able to retrieve 
the steady state solutions of the previous one-species model. The equivalence obtained between the effective carrying 
capacity and the extrinsic growth factor is complete only for a particular case, the Gompertz model. Here we unveil 
important aspects of sigmoid growth curves, which are relevant to growth processes and population dynamics. 
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1. Introduction 

Population dynamics models are useful when one 
tries to understand, describe or predict the behav- 
ior of a wide range of time dependent processes in 
several disciplines. To easily formulate and write 
the solutions of the population dynamics model, we 
present a one-parameter generalization of the loga- 
rithmic and exponential functions. This generaliza- 
tion has been first introduced in the context of non- 
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extensive statistical mechanics [1-5]. The general- 
ized logarithm or the so-called q-logarithm function 
is defined as: 

ln,(x)=. liin^^= /-^. (1) 

q'^q q' J t * 

1 

The natural logarithm function is retrieved for q — ^ 
0. The inverse of the ^-logarithm function is the q- 

exponential Junction 

{\mi{l + qxY'^' , \iqx> -1 
, (2) 
, otherwise 
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so that for q ~ 0, one retrieves the usual exponen- 
tial function. The use of these functions is conve- 
nient since it has allowed us to find and simplify (us- 
ing their properties) the solution of the models with 
time-dependent intrinsic and extrinsic growth rates. 

The simplest way to deal with population growth 
is to consider one-species models. In these models, 
individuals do not explicitly interact with the exter- 
nal ones and, at time t, the number of individuals 
is N{t) > 0, with initial condition A^o = ^(0) > 
0. The parameters are the intrinsic growth rate 
K > and the environment carrying capacity K = 
N{oo) > 0, which takes into account all possible 
interactions among individuals and resources [8] . If 
N{t)/K <^ 1, a general model is the von Foerster 
et al. model [9], where the per capita growth rate 
is dlnN/dt = kN^. Considering t = Kt > as 
a dimensionless time unity, the model solution is: 
N{t) = Noe-g[T/{qT)l with A^o = (gT)-i/9 and 
produces a population size divergence at a finite 
time T, obtained from the g-exponential qr/ (qT) > 
1 [10]. As g ^> 0, the exponential population growth 
[Malthus model) N{t) = iVge"*, with divergence at 
an infinite time, is obtained. The population size 
divergence can be dismantled considering a finite 
carrying capacity. 

The richness ecological community cannot be 
properly described only by the one-species models. 
The interactions between species becomes better 
formalized with the (Malthus-like) Lotka-Volterra 
equations [15], which only explains prey-predator 
behavior in its original formulation and presents 
stability problems. Besides the prcdation interac- 
tion, there are many other different kinds of inter- 
action taking place between two biological species. 
For instance, if species negatively affects each other, 
when they occupy the same ecological niche and use 
the same resources, there is competition. If species 
favor each other, there is mutualism or symbiosis. 
These ecological richness is partially appreciated 
in the competitive (Verhulst-like) Lotka-Volterra 
model studied in Rcf. [19]. This model presents 
complete analytical solutions. It has a non-trivial 
phase diagram and solves the stability problem of 
the Malthus-like two species model. 

Here, we consider a generalized two-species pop- 
ulation dynamic model and analytically solve it for 
particular ecological interactions. These two-species 
models can be simplified to a one-species model, 
with a time dependent extrinsic growth factor. With 
a one-species model with an effective carrying ca- 
pacity, one is able to retrieve the steady state solu- 



tions of the previous one-species model. The equiva- 
lence obtained between the effective carrying capac- 
ity and the extrinsic growth factor is complete only 
for the Gompcrtz model. 

Our presentation is organized as follows: In the 
Sec. 2, the Richards model is presented in terms 
of the generalized functions, which has the Ver- 
hulst and Gompcrtz models as particular cases. We 
show that the steady state solution of the Richards 
model with extrinsic growth rate is the same for 
the Richards model without extrinsic growth rate, 
but with a modified carrying capacity. Also, for a 
particular case of the Richards model — > 0), that 
is the Gompcrtz model, not only the steady state 
is the same but also the whole system evolution 
(transient). In the Sec. 3, we present a generalized 
model of two interacting species. We show that the 
interaction between species can also be interpreted 
as an extrinsic growth factor. This allows one to 
represent a two-species system by a one-species 
model with a modified carrying capacity. In Sec. 4, 
we present our final remarks. 

2. Generalized one-species model 

The growth of individual organisms [11], tu- 
mors [12] and other biological systems [13] are well 
described by mathematical models considering a 
finite carrying capacity which lead to sigmoid 
growth curves [14-16]. In this context, it is conve- 
nient to express the population size N with respect 
to its equilibrium value, i.e. p — N/K . The Richards 
model [17] 

-^lnp = ln^p, (3) 

binds the Gompertz (q ~ 0, dlnp/dr ~ Inp) and the 
Verhulst (g = 1, dlnp/dr = 1 — p) models through 
the q parameter, which microscopic interpretation is 
given in Refs. [4,18,21]. This parameter is related to 
the range of interaction between the individuals that 
compose the population and the fractal pattern of 
the population structure. This generalization is cor- 
roborated by several different approaches in terms 
of Tsallis statistics [6] and logistic equations of arbi- 
trary order [7] for the classical Verhulst and Malthus 
models. The solution of the Richards model is: 

P^^^l = ~~n — r// 1 -ri ' (^) 
eg{lnq[l/pQ]e ^} 

whose initial condition is po = ^(0) and the steady 
state solution is p* ~ p{oo) ~ 1 [20]. 
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Adding an extrinsic growth rate e to the Richards' 
model, one has the so caUed Richards -Schaefer 
model and respective solution: 



dlnp 
dT 

p{t) 



e,~{ln^[e,-(e)/po]e-[i+<^^l-} ' 



(5) 
(6) 



with steady state p(oo) = p* = eq{e), which behaves 
as an order parameter. The Richards' model solu- 
tion is retrieved for e = 0. The parameter e plays the 
role of an external factor that removes or insert indi- 
viduals in the population. For instance, it can be an 
interaction factor between different species, or the 
effect of the cancer treatment when one deals with 
tumor cells population. Survival of the species is ob- 
tained when qe > —1, and extinction occurs other- 
wise. The critical value e^'^-' = —1/q separates the 
extinction from the survival phase, which behaves as 



near the critical point. This tran- 



sition occurs for g 0, so that only the Gompertz 
model does not present the extinction-survival tran- 
sition. 

From the analytical solution of Eq. (6), we show 
that the insertion (e > 0) or removal (e < 0) of in- 
dividuals from the population modifies the steady 
state solution. Instead of considering the population 
size with respect to the medium (bare) carrying ca- 
pacity, we propose to consider the population size 
with respect to a rescaled carrying capacity K' = 
eq{e)K. It means that the N population individuals 
live in an environmental with an effective carrying 
capacity K' , but now without insertion or removal 
of individuals. The population size with respect to 
the new equilibrium value is p' ~ N/K' = p/eq{e) 
and one has the Richards' model and its respective 
solution: 



d\np' 

dT 

p'{t) 



-KM , 

p{t) _ 1 
eq{e) eg[lnq(l/p;,)e 



(7) 

(8) 



with steady state: p'(oo) — {p')* = 1, which of 
course leads to: p* = eq{e) in the original variable. 
Observe that the Eq. (7) is valid only for the sur- 
vival phase, i.e. qe > —1, otherwise K' vanishes. 

Comparing the solutions of Eqs. (6) and (8), one 
sees that the considered models present the same 
steady state solutions. Nevertheless, they have dif- 
ferent transient behavior, except for Gompertz 
model — > 0), where the equivalence between the 



original and rescaled models is complete. The plots 
of Fig. 1 illustrate this behavior. In fact, this prop- 
erty comes from the argument of the exponential 
function: while in Eq, (6) this argument is — t , in 
Eq. (8) the argument is —(1 -I- qe)T. Thus the evo- 
lution of the two systems become identical only if 
qe = 0. As we are dealing with e ^ 0, the evolution 
of the two model is the same only when g — > i.e. 
the Gompertz model. 




Fig. 1. Plot of population p(t) [Eqs. 6 and 8] as a function 
of T, with: po = O-lj f = —0.3, q = —1, g — > (Gompertz 
model) and q = 1 (Verhulst model). The full equivalence of 
the models with extrinsic growth factor and modified carry- 
ing capacity is only obtained for the steady state solutions. 
Nevertheless, for 5 — > 0, the transient solution is also equiv- 
alent. 



3. Generalized two-species model 

Here, we introduce a generalization of the 
Verhulst-like two species model, considered in 
Ref. [19], which we call the Richards-like two-species 
model. We consider the Richards term instead of 
the Malthus one in the original Lotka-Volterra 
equations. In this type of model, the two species 
interact according to the Lotka-Volterra equation, 
however the species also suffer from inter-species 
competition. As in Ref. [19], we let the interaction 
parameter vary from negative to positive values to 
cover all possible ecological regimes. The proposed 
model is written as: 



dpi 

dT 

dp2 

dT 



:pi[lngi pi -\-e1p2] 
--pP2[^nq^P2 + e2Pi] , 



(9) 
(10) 



where pi = Ni/Ki > 0, for i = 1, 2, where Ni > 0, 
Ki and Ki > arc the number of individuals (size). 
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net reproductive rate, and the carrying capacity of 
species i, respectively. The carrying capacity Ki rep- 
resents the restriction on resources that comes from 
any kind of external factors that do not have to do 
with species 2 and similarly for K2. Time is mea- 
sured with respect to the net reproductive rate of 
species 1, t = Kit > 0. The scaled time is positive 
since we take the initial condition at to = 0. More- 
over, the two net reproductive rates form a single pa- 
rameter p = K2/K1 > 0, fixing a second time scale: 
t' = pT ~ K-zt. In the right side of Eq. (9), the term 
Pi lugj^ pi represents the competition between indi- 
viduals of the same species (intraspccific competi- 
tion) and eipiP2 represents the interaction between 
individuals of different species (interspecific interac- 
tion) [15,16]. The right side of Eq. (10) behaves sim- 
ilarly. The non-dimensional population interaction 
parameters are ei and £2, which are not restricted 
and may represent different ecological interactions. 
Contrary to p, which has no major relevance to this 
model (since we consider only p > 0), the product 
£162 plays an important role, so that £i£2 < means 
prcdation; £i£2 = means commensalism, amensal- 
ism, or neutralism; and £i£2 > means cither mutu- 
alism or competition. For gi = (72 = 1, one retrieves 
the model of Ref. [19]. 

Here, wc restricted ourselves to the particular case 
when £i£2 = 0. Thus, only one species is affected 
by the other or they do not interact at all. Closed 
analytical solutions can be found in this particular 
case. Consider that the individuals of species 1 is 
unaffected by species 2. In this case, species 1 is 
described by Eqs. 5 and 6, so that: 



dlnpi 
dT 

Pi{t) 



= - liigi Pi 



.[ln9i(Pi,o)e 



(11) 
(12) 



where the initial condition is pi^o = Pi (0) and steady 
state: pi(oo) = pi = 1. The individuals of species 2 
are affected by species 1 through an extrinsic growth 
factor £2(7"), and one writes the equation for species 
2 as: 



d\np2 

dT 
where. 



P[ln,^2 P2 - ^(t)] 



<^2{t) = £2Pi(t) 



£2 



(13) 



(14) 



The solution of Eq. (13) has been obtained in 
Ref. [20]: 



P2{t) 



692 [4 (t)] 



fi [ e,-,[4(0)] \ e,- [4(0)1 



[l+92e2(-r)]pT 



where the initial condition isp2,o = P2(0) and steady 
state: ^2(00) = p^ = eq^{t2)- In Eq. (15), £2(t) = 
(l/r) /J" c?t'£2(t') is the mean value of £2(t) up to r. 

For £2 < (£2 > 0), one has amensalism (com- 
mensalism), so that species 1 adversely (positively) 
affects species 2. The bread mold penicillium is a 
common example of amensalism. The analytical 
solution for the amensalism and commensalism 
regimes is obtained considering the time dependent 
extrinsic growth factor as a function of species 1, ac- 
cording to Eq. (14). Penicillium secretes penicillin, 
a chemical that kills bacteria. Consider the two- 
species model of Eqs. (11) and (13), with £2 < 0, 
the interaction term e'2 can be interpreted as an 
extrinsic growth rate of a one-species model. As we 
have seen, an extrinsic factor can be incorporated 
into the carrying capacity. In this way, a two species 
model can be interpretated as a one-species model, 
with modified carrying capacity, as depicted in Fig 
2. 



Bacteria 
Carrying capacity - K 



Bacteria 
Carrying capacity - K' 




(a) 



(b) 



Fig. 2. Representation of the amensalism ceologieal inter- 
action between the bread mold and bacteria. A two-species 
model can be simplified to an one species model with a mod- 
ified carrying capacity (gray area), (a) Bacteria occupying 
all the available space with carrying capacity K. (b) Re- 
stricted bacteria growth due to the interaction with bread 
mold, carrying capacity K' < K. 

Species 2 evolves as a time dependent one-species 
(Richards-Schaefer) model, whose solution has 
been obtained in Ref. [20] and extinction occurs for 
£2(r)g2 < —1- The evolution of both species is pre- 
sented in Fig. (3) for two different values of qi. The 
time transient depends on species 1, through the 
effective extrinsic factor £2(t). For q^Q, there is a 
dependence on the mean value £2('''), see Eq. (15). 
Considering the Gompertz model (52 ^ 0), the 
solution can be simplified to: 



P2{t) 



oe2('^) 



P2,0 



4(0) 



j4(o)- 



,(^)-p^i 



(16) 
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5 10 15 
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(b) 

Fig. 3. For several values of q2, plots of pi(t) (thick lines), 
given by Eq. 12 and of P2(''") (thin lines), given by Eq. 15, 
where the integral of Eq. 15 has been calculated numerically. 
In both cases p = 1, ei = and £2 = 1/2: (a) gi = 1/2 and 
(b) qi = 2. Notice that if £2132 < -1, Pj = ^■ 

4. Conclusion 

We have proposed a generalized two-species pop- 
ulation model and have found its analytical solu- 
tion for the commensalism and amensalism regimes. 
From this result, we are able to show that the so- 
lution of the two-species model is equivalent when 
one considers one-species model with a time de- 
pendent extrinsic growth factor. Also, the extrinsic 
growth factor can be incorporated into the carrying 
capacity and the full equivalence, i.e. steady state 
and transient solution, is retrieved for the Gompertz 
model. Otherwise, the equivalence persists only for 
the steady state solutions. This feature reveals an 
important aspect of sigmoid growth curves. Consid- 
ering two well known growth models, Verhulst and 
Gompertz , besides presenting a similar sigmoid be- 



havior, they are differently affected by an extrinsic 
growth factor. For the Gompertz model, the extrin- 
sic growth factor can be fully incorporated into the 
carrying capacity; i.e. removing individuals from the 
population has the same effect if one considers that 
this population grows in a more limited environ- 
ment. Any kind of external influence can be seen as 
a modification in the environment limitations. Al- 
though the extrinsic factor can also be incorporated 
in the Verhulst model, it correctly describe only the 
steady state solution, the equivalence is lost when 
one considers the transient solution. 
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